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ABSTRACT 


This thesis is an examination of using a fielded doppler shift radar, specifically the 
AN TPQ-37 (FIREFINDER). to gather wind signature data for use in ballistic artillery 
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current rawinsonde data as thev affect the artillery probabilitv of kill P(K) against a 
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I. INTRODUCTION 


Meteorology data is used by all branches of the service to aid in performance of 
some portion of their operational nussion. This information 1s used by the artillery and 
air defense artillery for ballistic calculations. It 1s used by chenucal personnel to 
provide accurate down wind predictions for chemical and nuclear contaminants. The 
aviation industry has been searching for wavs to more accurately depict the current 
wind profiles (wind speed and direction at various altitudes) around major airports to 
reduce the number of severe wind related accidents. The National Oceanic and 
Atmospheric Adnunistration (NOAA) has conducted numerous tests and extensive 
empcummentation {Refs. |.2: pp. 48.25], where they conclude that radar can be used to 
measure the wind profiles in certain atmospheric conditions. NOAA has also 
eanemMecestiatsthesl.S. Amy: soA NYP O-37 (EFIREFINDER), an X-band radar, can 
provide wind profiles in the presence of precipitation, and at cloud heights, but can not 
furnish wind measurements routinely up to $ - 10 km in the optically clear atmosphere 
Ghemi@nll-msyo1. werclenoth {Kel 5: p. 3], Chapter | begins with an Executive 
Sumimary. which outlines the problem, the solution approach, the analysis, and 


conclusions. Chapter I concludes with an overall outline of the thesis. 


x) CUENECUTRE SUNINIARY 
1. Background 
a. Current Concerns 

Having accurate and timely meteorological data is verv important today 1f 
the probability of a first round hit on a target of known location by field artillery 1s to 
be iniproved. There are many factors that affect the accuracy of an artillerv weapon. 
The age of the tube (demonstrated bv the internal wear and overall effects of metal 
fatigue), the human errors that occur when aiming the tube (both in azimuth and 
elevation), and survey errors which occur when establishing the location of the gun are 
just a few of the non-meteorologically induced errors which affect the point of impact 
of the rounds. 

In the field artillerv, ballistic meteorology is concerned with deternuning 
atmospheric conditions in the area where artillery rounds or rockets will be fired. 


Atmospheric conditions along the trajectorv of a projectile directly affect its accuracy 


and may cause it to nuss the desired point of impact. The meteorologically imposed 
error can be as much as 5 to 10 percent of the range, even wmiderwsmioles earner 
conditions [Ref. 4: p. 2-1]. 

In the U.S. Army. meteorological (MET) data is used when computing the 
trajectories of artillery and mortar rounds, and when adjusting friendly fires. MET 
data is also used by sound ranging platoons in determining enemy artillery locations 
and by the chemical personnel when predicting the effects of the atmosphere on 
chenucal radiological contanunants that mav be introduced in the course of battle. 

For all of these uses, it is@mperatimemtial themeata, be as cunmentsamcmcas 
accurate as 1s physically and technically possible. 

b. Current System 

Currently. the U.S. Army relies on the Vieteoroiogicar Seciiommassiamcumrs 
the division artillery to gather. analyze, and distribute the necessary data. To do this, 
the MET section periodically flies weather balloons. The balloons carry a radiosonde 
instrument package which directly measures air temperature. humidity and barometric 
pressure. This data is continuously transmitted to a ground station which tracks the 
balloon as it rises through the atmosphere. The tracking data thereby provides wind 
speed and direction at the various altitudes (or zones) necessarv for ballhistic 
calculations using standard trigonometric and analytic geometry techniques. This data 
is analvzed and collated bv the MET section and then sent via secure radio-teletvpe to 
an operator located at the artillery battalion fire direction center (FO @yw mere 1tus 
entered into the tactical fire direction svstem (TACFIRE) computer bv hand. 
Henceforth. this data will be referred to as balloon MET data. This data is then used 
for all artillery firing calculations until the next meteorological data update. By 
doctrine, this 1s supposed to occur every four hours; however, in combat, where units 
are constantly on the move, it is verv doubtful that the update will occur that often. A 
more realistic figure nught be every six hours. This allows for the movement of the 
MET station, equipment setup, balloon preparation and flight, data analysis, and 
transnussion of the data to the FDC. 

This large time lag in data updates could potentially cause large errors in 
the accuracy of artillery fire nussions. Due to the dynamics of the atmosphere, this can 
happen in as little as two hours [Ref. 5], and in less time if a storm front moves 
through the area during that time, since more error is introduced into the data in these 


situations. 
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Cop piad: LET 

The requirement then is to update the MET data in a timely manner to 
ensure reasonable accuracy in artillery fire missions. One alternative is to accept the 
error and make the best of the svstem that is currently fielded. This should be 
considered as a baseline situation from which any improvement can be measured. 

Another alternative would be to double or triple the number of MET 
sections, thereby allowing for staggered balloon flights. accomplishing an update every 
two hours or so. The expense in both manpower and equipment quickly gets verv 
large with this alternative. 

‘A third alternative involves using a different svstem for obtaining MET 
data. One such system is a doppler shift radar with a phased array antenna which can 
provide the necessary wind profile when used in conjunction with a mathematical 
aeorithim, }Ref. 6], developed at the U.S. Armis Atinospheric Sciences Laboratory 
(ASL) for predicting temperature and barometric pressure at the various altitudes of 
the wind profile based on values obtained on the ground. The latter is known as an 
“analytic atmosphere.” The radar can be used as follows: at each altitude, the radar 
Sai eSmleralinesphereimeto directions, see Figure 1.1. The angle between the radar 
sampling beams should be about 70 - 90 degrees. Lower separation angles have been 
used but they did not appear to produce as accurate a wind profile [Ref. 3: p. 17]. 
After all altitudes have been sampled, signal processing techniques are emploved to 
extract the wind profile. Currently this signal processing is performed off-line from the 
radar, but software to accomplish real-time processing of the wind profile is being 
developed through efforts at ASL. 

They ialewatmiespiere is based on the current time of dav. current 
surface readings of teniperature and barometric pressure, and the munimum and 
niaXinium temperature readings from the previous day. It provides temperature and 
barometric pressure values at twenty-six altitudes which are then stored in a computer 
file along with the wind profile provided by the FIREFINDER radar to produce a 
computer MET message in a special format for input to the TACFIRE computer. 
Currently, this data is entered by hand, but because of the nature of the radar and the 
existence of communications links between the TACFIRE computer and the radar, the 
potential for automated MET data updating from the radar to the TACFIRE 
computer exists. The project at ASL to incorporate wind profiling into FIREFINDER 


capabilities along with the analytic atmosphere is known as WINDFINDER. 
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Figure 1.1 Radar Wind Profiling. 
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Or eas cOnemmeietiat the FIREPINDER radar can sense wind 
profiles on a limited basis. Whether the MET data provided by WINDFINDER is 
accurate enough for ballistic artillery has not been deternuned. 

2. Objective 

In late 1987, the US Army will start fielding an tmproved meteorological data 
Satnerimesss Stenmanowi as tic ie Mata Systent, (YDS), AN TMQ-31. This system 
1s reported to be a vastly improved radiosonde system, that tncorporates automatic 
data transmission to TACFIRE. The current fielding plan indicates that two MDS will 
be frelded at each of the infantry and armor divistons in the US Army. Field artillery 
brigades, and separate infantry and armor brigades are to receive one VIDS each. The 
Pee anemicide iillery Seiool] at Ft. Sill, OK will receive eight svstems for testing 
and training. and other systems will go into war reserves. In all, fiftv-ftve MIDS are to 
be procured at an estimated cost of $1.5 million each. An additional seven downsized 
systems are to be procured for the light infantry divisions. the 82nd ABN Division, and 
the 10lst ABN (Atr Assault) Drviston. 

This thesis is a pilot study whose objective is to determine, using the limited 
data available, if the FIREFINDER radar, operating in the WINDFINDER mode, 
can successfully augment the current balloon MET to provide data to the TACFIRE 
computer for ballistic artillery frre missions as a low cost alternative to the almost $100 
million \{DS svstem upgrade. 

5. Data 

Ideally, the way to compare two MET systems 1s to fire rounds under a wide 
range of conditions using each system (radar & balloon) and compare the accuracy of 
the results. Unfortunately, because of limited funds no actual firings have been 
condweted. 

During the period 28 Februarv 1986 - 8 March 1986, ASL gathered MET data 
at Yuma Proving Grounds, AZ. The data consisted of balloon MET and concurrent 
radar MET. The MET station at Yuma normally flies a balloon at O500hrs, OSOOhrs, 
1000hrs, 1200hrs and 1400hrs, unless funded projects have coordinated a change in that 
schedule. Due to a shortage of funds, the decision was made at ASL to gather radar 
MET data in conjunction with just the normal balloon schedule at Yuma. The results 
of the experiment were 19 balloon MET and radar MET data sets. Appendix A 
contains the first balloon MET and radar MET data sets in the format of a computer 


MET message. 
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4. Problem 

The only analytical tool readily available to perform ballistic artillery 
calculations as a function of meteorology input [Ref. 7], was a computer model 
developed at the US Army Ballistic Research Laboratorv (BRL), Aberdeen Proving 
Grounds, Marvland. One of the weaknesses of this model for the purposes of this 
thesis Was its deterministic nature. If a target was specified, then the model would 
provide a quadrant elevation and initial deflection for the gun to “hit the target.” 
Conversely, a quadrant elevation and initial deflectron could be specified and the model 
would indicate the location of the target at which the gun was “aiming.” So. using the 
niodel, there was no, direct way (o,¢onipare thestwo VIEW miata isets loc er iene 
model were used to determine the quadrant elevation and initial deflection using one 
MET data set (assunung that the aim point was the actual target), then switching 
Operating modes, the other MET data set could be used to determine a second “aim 
point” based on the previously determined quadrant elevation and deflection. The 
distance between this arm point and the target location would then provide a measure 
of the relative accuracy of the second data set to the firstg Phe model .as then misea 
with the balloon MET data set as the reference data set. The next problem: was to 
Geternune the absolute accuracy of the radar data using only the relative-to-balloon 
data that was generated using the BRL model. 

5. Model 

FM 6-I41-1 [Ref. 8: p. 5-6], specifies the probability of hitting, P(h), a target 
20 meters by 20 meters for various ranges. Using those values, it was a simple miatter 
to deternune the absolute accuracy of the balloon MET data that was necessary to 
produce the listed P(h) for each range. 

[laving the absolute accuracy for the balloon MET data, and the relative-to- 
balloon accuracy for the radar M[ET data, a model was required that described the 
interaction of those two parameters that would vield an estimate of the absolute 
accuracy of the radar MET data. 

If we assume that the target location (T) 1s known, then we can compare the 
two MET systems by comparing the results of using each system. Figure 1.2 shows the 
relative positions of the target and the two aim points. If the balloon MET 1s used, an 
alm point 1s produced. Ideally, this aim point should be located at T. [Iowever, the 
balloon MET system is not 100% accurate, so a random bias exists from T to the 


baltoon MET aim point B. This also applies to the radar MET system and its aim 


point R. The difference between T and B ts a random vector. Call it TB. Likewise, 
fom Pte Rois a random vector PR. The differciice of these two random vectors TB 
TR is another random vector RB. RB is the relative bias between the balloon MET 
afd the radar Weel and is included in the data generated by the BRL model, by virtue 


of using the balloon MET as a reference standard. 





Target 
i 
Balloon 
Radar B Am 
Aim R Point 
Point — 
y ite 
or 
/ r 
a 
i, 2 
/ ee 
te 
oe 


Figure 1.2 Target Location Model. 


6. MIOE 
In any analytical effort where a decision must be made between two or more 
alternatives, some measure of effectiveness (N[OE) must be used to keep the choice 
from being subjective. emotionally motivated and unsupportable in the face of test data 
and evaluations. Rockower [Ref. 9: p. 3], specifies that “a crucial part of the initial 
analysis is selection of the appropriate measure of effectiveness.” This MOE facilitates 
an objective, normally numerically valued, comparison between the alternatives. Since 


this was an analysis of the effect produced by two different MET data systems on 


ballistic artillery, the probability of kill. P(k). produced bx each Susteren. as.) eagerness 
Thee 
7. Analysis 

To represent the fact that targets are found at different ranges from the gun, 
three different ranges Were used in the analvsis. To eliminate any regional effects from 
prevailing winds. the P(k) was calculated for eight different directions of fire for each of 
the three ranges. These P(k)’s were then averaged to produce an average P(k) for the 
range distribution chosen. Based upon conversations with numerous US Army and US 
Marine Corps artillery officers at Ft. Ord@eCA. Fre. Sill. “OINeand= the wa 
Postgraduate School. Monterey, CA, three ranges, 10KmieiZkni amc ieee 
chosen as “realistic” ranges for which a 155mm howitzer would be emploved. Two 
range distributions were considered. one where the target was assumed to be equally 
likely at anv one of the three ranges, and a second distribution where the target was 
assumed to be at 12km half of the time, and a quarter of the time at 10km or 14km. 
These range distributions are purely subjective, and do not reflect any “expected.” or 
doctrinal distributions. 

The objective was to determine if the radar MET could augment the balloon 
MET, not replace it. Because the much lower cost of the radar’s software development 
compared with the alternative of procuring additional balloon MDS units was 
sufficient reason to choose the radar, so long as using the radar resulted in 


performance no worse than the alternative, the hvpothesis to be tested was: 
Hy: Radar MET is at least as good as stale balloon MET 

with an alternative of: 
Fj: Radar MET is not as good as stale balloon MET 

This translates into numerical symbologv as: 


Hy: Radar MET P(k) 2 stale balloon MET P(k) 
fH ): Radar MET P(k) < stale balloon MET P(k) 


There were 44 comparisons made for the analysis, see Table 1. Because of the 
linuted data, many of the comparisons have only a few data points, as indicated. All 
comparisons were made using the current balloon MET as the reference data set to 


provide quadrant elevation and deflection. 


| TABLE 1 
MET DATA COMPARISONS 


Comparisons Data Points 


Current radar 19 
“hr old balloon 
hr old balloon 
ir old balloon 
r old balloon 


hort vtoOnGsO 


; 
= 
Gj 
y 


I 
I 
‘| 
I 





[t was found that in all comparisons, with a level of significance of @ = .05, 

the alternative hypothesis (11,) could not be accepted. 
S. Conclusions and Recommendations 

Stnce the original data was so sparse and not gathered with this analvsis in 
und, results are not as conclusive as a properly designed experiment could show. 
Based upon the results, 1t can not be concluded that the radar MET 1s worse than the 
stale balloon MET. The large level of significance required to reject the hypothesis 
that the radar MET was as good as the balloon MET indicates the P(k) provided by 
the radar MET was statistically as good as that provided by the balloon MET. As 
such, the conclusion that the radar MET can be used to augment the balloon system 1s 
very strongly supported. 

Given such strong indications of parity between the two MET systems, the 
next phase of testing should be pursued. These results warrant the expense of a field 
test Where both MET svstems are employed and targets are engaged with artillery of all 
calibers. The test should employ the balloon tn its doctrinal operational mode, and the 
radar should be emploved every two hours. The guns should fire under the same target 
conditions and should employ the MET data available at that time. The comparison 


in this case must be of the actual miss distance from the target to the point of impact. 


5B. SHESIS OVILINE 
The outline for the remainder of this thesis is basically the same as that used for 


the Executive Summary. The areas summarized above are discussed in much greater 


depth in the respective chapters. In Chapter II, the data gathering and preparation is 
described in more detail. Following that. Chapter III presents a development of the 
analvuc model. Chapter IV discusses the selection and calculation of the MOE, while 
the calculations necessary to deternune the P(k) are developed in Chapter V. Chapter 
VI presents the analysis of the results of the tvomnamce distributions Gates 0 
contains the final conclusions and recommendations. Appendix A contains an example 
of the balloon MET and radar MET data sets. Appendix B has values for the results 


of each individual range. 
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If. DATA PREPARATION 


The original \fET data was collected at Yuma Proving Grounds, AZ by ASL, as 
described in Chapter [. Due to the sparseness of the data. a method of comparing the 
effects of the radar MET data against the balloon MET data had to be devised. The 
only analytical tool to compare MET effects on artillery fires which was readilv 
available and accepted throughout the artillery and ballistic meteorology communities 
(Ref 7} was the General {rajectory Prosram (GIRAJ) written by the US Army 


Ballistic Research Laboratory, Aberdeen Proving Grounds. MD. 


A. GTRAJ 

The General Trajectory Program was developed to assist in the development of 
artillery animunition, propellants. and related products [Ref. 10]. It is written in 
Benlivy..// and ainstalicd on a DEC PDP Ii/SO (VAN) series computer. It uses 
the point mass or modified point mass trajectory miodel to compute trajectories for 
boosted and non-boosted projectiles [Ref. 11]. As such, the model is completely 
Geternunistic and does not incorporate ballistic dispersion in its output [Ref. 12]. 
Washburn [Ref. 13: p. 1], describes ballistic dispersion as a probability density function 
for firing errors which incorporates the round to round variation in both cross-range 
and down-range impact point values. Since no ballistic dispersion is assumed, the 
output from GTRAJ is where the gun was “aimed.” given the initial conditions 
assunied and the MET data base in use at the time [Ref. II]. This deternunes the aim 
bias. [t is this aim point that will be used to calculate a P(k) after incorporating 
ballistic dispersion into the analysis. 

GTRAJ is a fully menu driven program that can be operated either interactively 
or from a batch file. It has many different integration modes that may be specified 
[Ref. 11]. The two used for this analysis are 1) integrate backward from a target of 
Known location to determine the necessary quadrant elevation and deflection to aim 
directly at the target, and 2) given a quadrant elevation and deflection, integrate 
forward to determine the perceived aim point. Using this model, there was no way to 
direct) evallate the independent effects of both MIET data sets. Howeveryaf the 
model was used with one MET data set to determune the quadrant elevation and inital 


deflection necessarv to aim at a known target, then switched to the other operating 


mode. the other MET data set could be used to deternune a perceived aim point using 
the previously determined quadrant elevation and deflection. This aim point would 
not, in most cases, be the same as the originally chosen target location. This aim point 
would then provide a relative accuracy of the second MET data set compared to the 
first. The current balloon MET data set was always used in conjunction with the first 
niode of operation to deternune the quadrant elevation and deflection necessary to aim 
the gun at the target. Using this output with the corresponding comparison MET data 


set, the second mode was used to determine the corresponding aim point. 


B. GENERATION PROCESS 

‘A copy of GIRAJ was provided by ASL and installed on a DEC PDP I1 780 
(VAX) located in the Computer Science Department at thle “Naval Posteraduare selaer 
Wioiiterey se 

The generation process 1s depicted at Figure 2.1. The model must be initialized 
by specifving the weapon system, amimunition type, mode of input (interactive or 
batch}. and output file name. Next, the mitial conditions are selected, the azimuth is 
set to 0.0 mils, and the mode of operation 1s set to integrate backward to determine the 
quadrant elevation and deflection for the given range. The balloon MET data file is 
specified and the computational parameters are set. This initial output produced is a 
quadrant elevation and deflection. Next, the operating mode 1s changed to integrate 
forward, the radar MET data file is specified, and the quadrant elevation and deflection 
just calculated are input. Keeping all other input parameters the same, GIRAJ is run 
again. This output is a range and deflection which is the perceived aim point for the 
radar MET in relation to the aim point of the balloon MET (target location). This 
point is stored in a separate data file for further processing. The next step is to 
increment the azimuth of fire by 800 mils (45°). The balloon MET data file is again 
specified. operational mode | (integrate backward) is selected and the process starts 
over again to generate another radar MET aim point. This process continues until aim 
points have been generated for the eight directions (N,NE.E,etc) for each of the three 
ranges, 1Okm, 12km, and I4km. All comparisons were made using the current balloon 
MET as the reference data set to provide quadrant elevation and deflection. The 
inputs that are required by GIRAJ [Ref. 11], and the values used for this analysis are 


in Table 2. 
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Figure 2.1. Data Generation Flow-chart. 
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Required Input 

Initialization 
Weapon Caliber 
Weapon system 
Ammunition 

Initial conditions 
Azimuth of fire 
Curved earth or flat 
Locationyof guns. 2) 
Location of target 
Range to target 


Mode of operation 


Computational parameters 
Charge 
Muzzle Velocity 
Powder Temperature 
Quadrant Elevation 
Deflection 
Weight of round 
MET data file 


Value Used 


155mm 
MI09A1 
NOG eel le 


Variable 

Curved 

O00 

Variable, Variable,O 
10km. 12km. 14km 


Integrate from 
target to gun 


Integrate from 
gun to target 


TV 
568m $ 
70° F 
Variable 
Variable 
95lbs 
Variable 


C. REFERENCE ORIGIN 

The output from GIRAJ is in the form of a4 ranee@itom the oun tortie poinucr 
aim and a deflection, (D), which is a distance left (-) or right (+) of, and perpendicular 
to, the gun — target line. To know the exact point of aim, it is necessary to also know 


the direction of fire. To facilitate automated calculation of P(k), this aim point was 


converted to an \.\ coordinate using the target location as the origin. In this manner. 
the dependence upon the direction of fire was removed. A standard translation of axis, 
as mentioned tn Thomas, [Ref. I4: p. 518], was used to transform this orrginal arm 
point, which was a function of range to the target and direction of fire. to an N.Y 
coordinate independent of the directron of fre. The conversion to the reference origin 
ii oOlwes (Ntcempan-es wm cictringsto Figure 2.2, locate the target m the firing plane 
Where the gun 1s assumed to be at the orrgin. The target ts located at the pornt 
(X).¥ 1). This is a standard polar to rectangular coordinate conversion process, where 
R, is the range to the target from the gun. and § ts the direction of fire tn degrees, 
using the standard rectangular coordmate svstem where the abscissa is associated with 
Us 


A= R,eCos@ 
2] = R,eSin@ 


Next. the aim point must be located in this plane relative to the origin. The point 
(X5.¥ 4) 1s the location of the aim point. where R95 is the distance from the gun to the 
ain. potnt, and D, the deflection. is the perpendicular distance left (-) or right (+) of 


the gun — target line. 


V(R>2— D*)eCos8 + DeSind 


a 
2 
ll 


Y, = ¥(R>?—D2)eSinO — DeCos6 


Lastly, the point on) must be converted to a point (X,Y) relative to Ge) Which 


will be the new origin. 


Y= ieee Y | 


Once the data has been converted to this reference origin, it is much easier to automate 
the calculation of the P(K) for that range and MET comparison. 

Prior to generating the data, it was necessarv to develop a model that would 
explain the relation of the radar aim: points with the balloon atm points. This ts 


discussed in the next chapter. 


CAG aa 


Sar) 





Ny = R,eCos8 

Yy =~ Ry e¢Sin8 

X,= (R25 —D2)eCos® + DeSin® 
Yo = ¥(R2, — D~)eSinO — DeCos8 
Y= Ne am Yy 


Figure 2.2. Transformation to Reference Origin. 
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Ll fRODEL DEVELOPMENT 


With anv problem. there is more than one approach which will culmunate in a 
solution. With this problem. the methodology of the analvsis 1s very important. Due 
to the complete lack of live firings. the linuted amount of initial MET data available. 
and the fact that only this MET data was used to produce the data for the analvsis, it 
Was important to choose an approach which would maxinuze the usefulness of the 
MET data. This chapter describes the model used in the analysis. The first section 
defines the problem that the model had to fit. The second section describes the 


deternunation of the different elements of the model. 


A. MODEL 

-As discussed in Chapter I]. the linuted amount of data and the method of 
emploving the computer model GIRAJ required an analvtic model that could describe 
the absolute accuracy of the radar MET when what was being measured was the 
accuracy of the radar MET. relative to the balloon MET. 

To construct the analvtical model. the target location (T) 1s assumed to be 
known. If balloon MET is used. an aim point for the artillery piece is determined. 
This aim point describes the random vector TB, see Figure 3.1. from the target to the 
aim: point. Ideally, this aim point should be located at T. However, MET data has 
random errors associated with it that are approximately constant from shot to shot. 
hence TB is a random variable. and B will not normally coincide with T. TB is 
distriouted about T with some mean. Hp. and variance, aa For consenience of 
calculation. TB is assumed to be distributed according to a circular normal (0.07) 
distribution. The mean is assumed to be Q. as any consistent error recognized during 
testing of the system, could be compensated for during calibration of the svstem. 
Likewise. if radar MET is used. TR. another random vector is generated, also 
distributed about T, such that TR ~ Organ) TR and TB are independent random 
variables. There is one additional random variable in the model. RB. RB 1s the 
difference between TR and TB. Since RB is the difference of two normally distributed 
random variables, it is also normally distributed, [Ref. 15: p. 267], RB ~ Gan): 
werrere oe = oe *- e In this analysis, it is the variances of these random 


variables that are of the most concern. 
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Figure 3.0 Target Location Model. 


B. BIAS DETERMINATION 
‘A measure of pp is implicitly contained in the output aim points of GIRAJ 


2 


— : ae ? 
as it 1s the radar-to-balloon bias. Consequently, if either 6°p or Op could be 


calculated or deduced, using the fact that Go? 


RB can be estimated with GIRAJ, the 
third variance could be discovered. 

Through many tests and experiments, the US Army has determined the 
probability of hitting, P(h), a target of two different sizes, lOm x 10m and 20m x 20m. 
These probabilities are specified in FM 6-141-1, [Ref. 8: p. 5-6], and reproduced in part 


in Table 3. 


26 


The probability of hitting a target uniformly distributed across an area A, 1s 


given by Equation 3.1. The firing distribution is f{x.v). 


P(h) =f f f{x.y)dxdy (eqn 3.1) 
A 


If we assume that X 1s independent of Y. f{{x.y) can be represented as f(x) @ f(y). Then 
E@uatOm.. 1 Pecones E@liation 5.2, where the target 1s assumed to be centered at the 


origin and has dimensions of 2a x 2b. 


a b 
P(h) = \ f(x)dx@ { ody (ean. 3-2) 
-a -b 


The US Army makes the assumption of bivariate normal firing errors [Ref. 16: p. 
13-1]. and has published the values of cross-range and down-range ballistic dispersion, 
Gy and Gy, that have been determined empirically. [Ref. 17: p. A-5]. Table 4 contains 
an extract of these values. By setting P(h) equal to the published value for each range, 
Equation 3.3 and a set of tables for the standard normal distribution, as found in 


[Ref. 15: p. 580], can be used to iteratively solve for Op. Column 4 in Table 3 is the 


published P(h), column 3 is the calculated P(h) given the inferred value of Op listed in 





colunin 2. 
x2 Cm ce 
P(h) = expl- dx # ——— | exp/- ay eqn 3.3) 
(h) ae Sa a) oo > Jdy — (eqn 3.3 


where Ll = 2no7 | lee 2no-, 37, = g*y rs 67, , and a”, = oy - op. 


Once G75 has been determined, cane can be calculated. Knowing G5 and GRR. 


Pies eo: 
G~p is then determined by Equation 3.4. 


2p = oath — 57, (eqn 3.4) 





TA GIRS 
BALLOON BIAS VS RANGE T@sh Ge 


Target Size 
10min 
RANGE  BIAS(¢p) P(h) PM 6-I41-1 P(h) 


1Okm 37 505 OS] 
12km 18.16 03395 O34 
I4km 20.38 26s Ze 


Target Size 


20m x 20m 
RANGE  BIAS(op) P(h) FM 6-141-1 P(h) 
1Okm 16.09 Ose .166 
ie 20.01 mleOS alee 


l4km de-Ol 0928 





TABLE 4 
BALLISTIC DISREKSt Om Be nES FOR 155. Ow Ze: 


RANGE Oy Oy 
10km 8 a9 
12km 10 44 
14km 13 S1 


Mathematically, it makes no difference whether the arbitrary origin is located at 


the target. or if it is at the aim point provided by the balloon MET. Because it is more 


convenient in using GIRAJ, the remaining analvsis will assume that the balloon MET 
aim point is the origin, and BT 1s distributed about B; BT ~ \(0,67 8). 

After selecting the appropriate MOE, which will be discussed in the next chapter, 
the calculated values of op and cae will be used to generate the data to be analyzed 


in Chapter VI. 


IV. MOE DETERMINATION 


Chuvev [Ref. 18: p. 11], emphasizes that any MOE must be representative of the 
decision required, be simple in application, and contain, if possible, all of the basic 
elements under study. [n this thesis, two different methods of meteorology data 
collection are compared. The element of interest is how they individually affect the 
accuracy of ballistic field artillerv, as it impacts on mission accomplishment, when the 
artillerv uses each of the atmospheric models for targeting calculations. As such, the 


VMIOE should reflect that interest. 


A. ALTERNATIVES 

Ideally, to accurately compare the two MET systems, an experiment should be 
conducted, where numerous artillery fire mussions of all calibers are fired. Such an 
experiment would be very costly in terms of equipment, ammunition expenditures, and 
man hours involved in the conduct of the experiment. In such an experiment, the miss 
distance from the target would be an integral part of the MOE. In this analvsis, the 
distance from the target to the aim point could be considered as a “muss distance.” 
This would be a viable alternative if the firing distribution were circular normal. We 
want to choose an MOE closely identified with mission accomplishment of field 
artillery, 1.e. destroving the target. Hence, the probability of kill, P(k). incorporates 
artillery “accuracy,” as a function of the MET data, (represented by Op or Op), in that 
non-linear formulation most closely associated with mission objective. Hence, the 
higher the P(k), the better the aim pomts) [he approach thei isstosealetmatentiicat (ie 
produced by each system. This approach has intuitive appeal, and it has a closed form 


solution under the proper assumptions. 


B. CALCULATION OF MOE 

Once the MOE has been chosen, it becomes important to determine how to 
calculate the MOE. Since the radar MET is extrapolated above the planetary 
boundary laver (PBL), and balloon MET is measured [Ref. 5}, the problem 1s to 
determine how to remove the effect (bias) of the prevailing winds above the PBL. The 
method chosen was to “fire” at a target of given range from eight different directions, 
(N.NE,E.SE,etc.). Therefore, this firing method would have the effect of mitigating the 


bias if the maximum) ordinate of the round was above the PBL. 


30 


Assuming each of the eight directions of fire are equally likelv, the calculation of 
the MOE proceeds as follows. Determune each aim point for the eight directions and 
calculate the P(k) for each ain: point. The P(k)’s were then averaged to deternune a 
point estimate for the P(k) for that range, Equation 4.1. This was repeated for each 


range of interest. 


I 
range og & 


P(k) PUK) direction (eqn 4.1) 


The comparison P(k), which is the average of the range P(k)’s could be 
deternuned in many different ways. The method used should reflect the suspected 
range distribution of the targets engaged. As such, it could be viewed as a weighted 
sum of the range P(k)’s, Equation 4.2, where the weights were assigned based upon the 
suspected distribution of the ranges. For this analysis, two discrete distributions over 
10km. 12km, and I4km were chosen. The weights were the probabilities of the range 


osenmime,. seevlable 5. 
P(R= > Wenge ® P(K) (eqn 4.2) 


‘range 


TABLE 5 
TARGET RANGE DISTRIBUTION 


RANGE DIST 1 (p) DIST 2 (p) 
10km 333 .250 
12km 333 500 
Idkm 333 .250 


The next chapter presents the detailed calculation of the weighted P(k), the MOE 
selected for the system comparisons. This appears to be the most logical MOE, and 1s 


readily calculable. 


V. CALCULATION OF ri) 


When calculating P(k), the items of concern are target location, where the 
weapon is aimed, the firing distribution, and the distribution of the lethal effect of the 
ordnance being fired at the target. For this analysis, the target was assumed to be a 


point target located at the origin in the XY plane. 


A. DAMAGE FUNCTION 

The probability of destroying a target is the product of the probability of hitting 
the target. P(h), and the conditional probability of killing the target given that it is hit, 
D(r). commonly referred to as the damage function. As explained in Eckler and Burr. 
[Ref. 19: p. 16], if we assume that a target is located at the origin (0,0) in the XY plane, 
we can denote the probability density function of weapon impact points bv f(x.v) and 
the probability of destroving the target if the weapon impacts at (x.y) by D(x.y). Then 
the unconditional probability of destroying the target with a single round 1s given bv 


Equation 5.1. 
P(k) = J{D(xy)ef(x.y)dxdy (eqn 5.1) 


The damage function, D(x.y) is actually a conditional kill probabilitv. Although 
not required, in general the damage function is assumed to possess circular symmetry 
and be non-increasing, [Ref. 13: p. 2]. This means that the damage function is a 
function of only one variable, r = (X24 7) he. Since it is a probability, it ranges 
from. one to zero as r increases from zero to some maximum lethal radius, R, away 
from the target. The damage function represents the probability that R 1s greater than 


r. As such, 
D(r) = P(R > r) (eqn 5.2) 


for this analysis, two damage functions were considered, the so-called cookie cutter 


weapon and the diffuse Gaussian or Carlton weapon. 


1. Cookie cutter 

Since the damage function is the conditional kill probability, the cookie cutter 
weapon can be easily visualized as a lethal raditts, R. around the target. This radius is 
constant and forms a circle around the target, hence the term cookie cutter. If the 
Weapon strikes within this constant radius. r S R, the target will be destroved. 
otherwise, the target 1s undamaged. Conceptually. this weapon has a great deal of 
appeal. However, as Washburn describes [Ref. 13: pp. 4]. when the firing errors are not 
circular normal. or the aim point is offset from the origin, there is no closed form 
solution for P(k), and numerical integration or other numerical techniques must be 
emploved. 

Both of the above situations existed in this analysis. the ain point is not 
located at the origin and the probability density function for the firing errors, although 
bivariate normal. was not circular. 

-. Diffuse Gaussian 

The diffuse Gaussian weapon is one of the alternative damage functions that 
allow a closed form solution under the above conditions. The diffuse Gaussian weapon 
does not assume that the lethal radius is a constant. Instead. r is a continuous random 
variable and has a range of (0,©). The damage function for the diffuse Gaussian has 


the form of Equation 5.3 for some scaling factor, b. 
Lee 5.3 

D(r) = exp(—r*, 2b*) (eqn 5.3) 
The lethal area for this weapon then becomes 2nb?. During this analysis, b was 
chosen such that the weapon lethal area was equal to the lethal area covered bv a 
135mm high explosive (HE) round that has a bursting radius of 50 meters. Therefore, 
Equation 5.4 describes the relationship between the damage function and the bursting 
radius of the 155mm HE round. 


2nb? = mR? (eqn 5.4) 


Setting R equal to 50m, and solving for b, yields b = 35.36. 


B. FIRING DISTRIBUTION 

The US Army has, over the vears, conducted numerous tests to deterniune the 
firing distribution and its parametric values for the indirect fire weapons in the arsenal. 
The distribution specified in DARCOM-—P 706-101, [Ref. 16: pp. 3-1.A-5], 1s an 
elliptical bivariate normal distribution where both the cross-range and down-range 
components of the ballistic dispersion vary with the range to the target. The firing 
distribution is shown in Equation 5.5, where (X,.Y,) is the aim point, and (X.Y) is the 


actual point of impact, with Oy and Oy the ballistic dispersion. 


» » 
l NX-X,)* oy). 
figura ree ae santa ae cael }) (eqn 5.5) 
LNG Oy 10° y 20" 


The values for Oy and oy as found in FM 101-61-5-3, [Ref. 17: p. A-5], are extracted 
and appear in Table 6. 


TABLE 6 
BAUCLIISTI@ DISPeIsTo noice FOR 155 viv HOWIZER, 





C. DETERMINATION OF P(K) 
1. Numerical Determination 
The assumptions of normal firing errors and the diffuse Gaussian damage 
function combine very nicely to produce a closed form solution for P(k). Using 
Equations 5.4 and 5.5 in Equation 5.1, we obtain Equation 5.6. As Washburn, 


[Ref. 13: p. 5], states, this equation assumes the center of the error distribution 1s 
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(Hay ply-) with ballistic dispersion of (6x ,6,-). For this situation, the target is assumed 


to be located at the origin with certainty. 


, b- Hy "7: Z 
ae 6 aa 5 Ga) Ke 5-6) 
V(bo+ Ose )(- + 6-7) eC) 2( bea oa) 

As described in Chapter II, the way GIRAJ was used required one of the 
MET data sets to be used as a reference during the comparisons. Since the balloon 
MIET was the reference, all of the comparison aim points assumed that the balloon 
MET aim point was the origin, and the target was a random distance away from the 
balloon MET aim point. This does not invalidate any of the equations developed so 
far. but the target bias. Op. must be incorporated into Equation 5.6 to produce the 
probability of killing the target. Thus Equation 5.7 is the probability of killing the 


target with a single round, where (X,.¥,) is the GIRAJ generated aim point. 


- 2 , 2 
b-exp a No a a \ 
2 b- +67. +67) b> +07,-+ cael | 


P(kK) = (eqn) 


2 Awr+er+02p)! 2 
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When applying Equation 5.7 to the aim point data, it was necessary to recall 
that in the case where the radar was being compared against the current balloon. ome 
was implicitly contained in the offset aim point. As such, to calculate the absolute 
P(k) for the radar MET, it was necessarv to subtract out the balloon bias when making 
the calculations. 

2. Graphical Determination 

Figure 5.1 illustrates a graphical method of determining the radar bias, Op, 
which was used to verify the calculations described above. 

The plots represent the change in the arithmetic mean of the P(k) for each 
type comparison, as the value of op increases from 0 to 100. The plot of the current 
balloon (top line) can be thought of as an ideal standard, as for that case, the aim 
point was the origin and the only errors were due to ballistic dispersion. This 1s 


geation >./wmith «, = ¥Y 0, and op varying along the horizontal aXIS as 


O 
indicated. 


30 


SS = a — 
a 


I & DIRECTION) 





a St 
nr nt 
Se 


1 
I 


PK(AVG OVER RAN( 


PK VS TARGET BIAS 


0.6 


Ee Gin ap) 
CURRENT BALLOON 
CURRENT Ses 


O23 


f 


Q. 


cs 


0 





0.0 


0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0 
BIAS(Op) 





Figure 5.1 Graphical Bias Determination. 


To determine the total amount of bias represented by the average radar MET 
P(K). extend a horizontal line from the radar \{ET curve at the ordinate to intersect 
Witae@me ideal curve. Then drop a vertical line to the abscissa. The point of 
intersection is the amount of bias in addition to the ballistic dispersion necessary to 
achieve that P(k). The solid line represents this process on the radar P(k). 

Since the value of Opp is represented implicitly in the aim point for the radar 
ne leit 1s mecessary to subtract the balloon bias Op so the P(k) will reflect the 
absolute radar accuracy. 

Micevatesor (mc savenage tadar VIET P(k)-on Figure 3.1 is .2164. This 
represents an average bias of approximately 60.5 meters as read from the graph. This 
bias is the radar to balloon bias Opp. If we substitute the average Op (19.57) 
represented in Table 3, into Equation 3.4, an average radar bias of approximately 57.25 
1§ calculated. To determine the absolute average radar P(k), enter the graph at the 
bottom with this number (dashed line) and reverse the process described above. This 
vields an average radar P(k) of approximately .23. 

The remaining four curves on the graph describe balloon MET P(k)’s at 
various ages. These curves were calculated assuming 6p = 0. As such. to reflect the 
“absolute” balloon MET accuracy, the value of 6p must be added as in Equation 5.7 to 
produce accurate balloon MET P(k)’s. Graphically, enter the graph at Op = 19.57, 
and go up until intersecting each balloon MET P(k) curve. Then extend a horizontal 
newer the Ordinate to determine the absolute balloon VIET P(k) for each curve. The 


relative and graphical absolute values of the P(k)’s are included in Table 7. 


ae tee 
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As is apparent from Figure 3.1, when the balloon MET bias was incorporated 
into the results, the radar MET P(k) got better since the balloon MET was subtracted. 
and the P(k) provided by the stale balloon MET got worse as the balloon MET bias 
was added. 

Since the individual radar aim points inherently contained the value of Opp. a 
way of extracting the effect of Op. thereby leaving only the effect of Op, was required. 
Recognizing that the target distribution was circular normal, the length of the vector 
produced from the origin to the aim point could be scaled to reflect the “magnitude” of 
the radar bias. To accomphsh this scaling, each aim point produced in the radar 


comparisons was transformed as in Equation 5.8 
as. 4 ) ? : 
OS Yo) = VC pace cin «eam (eqn 5.8) 


These new aim points were used in Equation 5.7 to calculate the individual point 
P(k)’s. These P(k)’s were averaged to produce a P(k) for the given range, and the 
range P(k)’s were combined as described in Chapter IV to produce the comparison P(k) 


used in the analysis in the next chapter. 
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VAN vi eSIssOr RESULTS 


This chapter is a detailed analysis of the results of the comparisons in Table 1. 
The analvsis is both graphical and statistical. Section | 1s concerned with the graphical 
analvsis, Section 2 follows up with testing of statistical hypothesis. In this chapter, 
only the data from the two range distributions 1s considered. The individual range data 


and some ancillary graphical displays are included in Appendix B. 


A. GRAPHICAL ANALYSIS 

The graphical analvsis consists of two parts, a review of the data histograms as 
compared to a normal density function with the results of the Kolmogorov-Snurnov, 
(K-S), goodness of fit test. [Ref. 20: p. 346]. for a normal distribution, and analysis of 
the empirical quantile - quantile (Q — Q) plots. 

As indicated earlier, the first range distribution assumed that a target had an 
equal probability of being located at either 10km, 12km, or 14km from the gun. The 
only comparisons considered in this portion of the analysis are those with five or more 
data points. 

As is evident from the histograms in Figures 6.1 — 6.3, the radar MET P(k) data 
is the only one which remotely resembles unimodality. None of the histograms suggest 
similarity of distributions. However, due to the small sample sizes, similarities could 
not be ruled out. As indicated in Table 8, each of the samples has a rather high level 
of significance for the K-S test statistic, which indicates a reasonably good fit when 
compared with a normal distribution, even with these few data points. The radar 
comparison was the only sample which had enough data points to allow a Chi-Square 
goodness of fit test to be performed [Ref. 20: p. 189]. 

As found in Chambers, et al., [Ref. 21: p. 68], the empirical Q — Q plot can be 
used to determine if the two data sets differ by an additive or multiplicative constant. 
The points plotted in Figures 6.4 and 6.5 are the values of P(k) associated with every 
Pie eciicle cee > ll), =-.,.95), Of each data set. Each set of quantiles was 
calculated based upon the individual sample, i.e. given the sample, the 5th, 10th, etc. 
quantile was determined. The quantiles for the two hour or four hour balloon MET 
were then plotted against the radar MET quantiles. As seen in Figure 6.4, the majority 


of the plotted values lie above the X=Y line (solid line). A line parallel to that line 
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(dashed line) can be drawn such that approximately half of the data points are on 
either side of it. The dashed line represents a shifting of the N=Y line by 
approximately .05 units. This indicates that on the average, the two hour MET P(k)’s 
are approximately .05 higher in value than the radar MET P(k)’s. In Figure 6.5 it is 
possible to construct a line shifted by approximately .02 units. Thus, the four hour 
P(k) have a value that 1s approximately .02 higher then the radar MET P(k). As ts 


shown in the next section, this is not necessarily statistically significant. 


B. STATISTICAL TESTING 

Although all of the data sets had high K-S significance levels for a normal 
distribution, two transforms suggested by Bartlett, (Ref. 22: p.52], were applied against 
them; the so called Fisher's Z transform where Z = 1/2 In{(1+p) (1-p)], and the arcsin 
transform where Z = arcsin(./p). Assuming that the goal was to spread the data and 
stabilize the variance, a change of variable was imposed where r = 1-2p, then the 
lisher’s Z transform was used on r. This yielded Z = 1/2 In{(l-p) p]. But since it was 
more appealing to have large values of p map into large values of Z, and likewise for 
small values of p, the reciprocal, Z = 2 In[p,(1-p)], was chosen. None of these 
transforms produced anv increase in the K-S significance levels. 

The objective was to determine if the radar MET could augment the balloon 


MET, not replace it. Because the much lower cost of developing, testing, and fielding 
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Figure 6.1 Histograms of Radar P(k). 
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Figure 6.2. Histograms of Two Hour Balloon P(k). 
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Figure 6.3. Histograms of Four Hour Balloon P(k). 
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Figure 6.4 Q- Q Plots Radar P(k) Quantiles vs Two Hour P(k) Quantiles. 
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Figure 6.5 Q-Q Plots Radar P(k) Quantiles vs Four Hour P(k) Quantiles. 





the radars software compared with the alternative of procuring additional balloon 
MDS units was suflicient reason to choose the radar, so long as using the radar 
resulted in performance no worse than the alternative, the hypothesis to be tested was 


chosen to be: 

Hy: Radar MET is at least as good as stale balloon MET 
with an alternative of 

H): Radar MET is not as good as stale balloon MET 
This translates into numerical svmbology as: 


Hp: Radar MET P(k) 2 stale balloon MET P(k) 
I}: Radar MET P(k) < stale balloon MET P(k) 


It was apparent that as the data sets got larger, the more closely they fit a 
normal distribution. Although the data sets “fit” a normal distribution, there was some 
hesitancy to apply the two sample t-test with impunity due to the small sample sizes. 
As such, in the interest of conservative estimates, both parametric and non-parametric 
tests Were performed. The data produced by the comparisons listed in Table 1 are 
shown in Table 9, the results for the individual ranges are in Appendix B. This data 
Was tested using the two sample t-test with Welch’s approximation to the Behrens- 
Fisher problem of unequal variances, [Ref. 15: p. 451], and the Mann-Whitnev test, 
[Ref. 20: p. 216] to determune if the data came from the same population. The Mann- 
Whitney non-parametric test was chosen because it does not assume that the data sets 
come from any underlying distribution, whereas the two sample t-test assumes that the 
underlving population is normally distributed. The results of those tests are 
summarized in Tables 10 and It. Both tests were performed using MINITAB, a 
commercially available computer statistical analysis package. 

As is evident from the tables, neither of the tests resulted in rejection of the null 
hypothesis and acceptance of the alternative hypothesis. As such it can not be 
concluded that the P(k) provided by the radar MET is any worse than the P(k) 
provided by the stale balloon MET. There was some hesitancy in testing the 5hr, Ohr, 
7hr, and 9hr data individually, since they each had only two data points, however, by 
combining them into one data set, a test for balloon data over four hours old was 
conducted. This is the last line in each distribution in the table. The value of a’ is the 


level of significance that would have to be used in order to reject the null hypothesis. 
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TABLE 9 
CNE@er rel) Fen) FOR Te? CONIPARISONS 


Ree\GE Oils PRIBL TION =1 


RADAR Divs Sa SHR 
0.128655 0.324017 0.402227 0.335361 
0.230003 0.451937 O.32941] 0.366009 
0.200316 O.137 741 0.240176 0.104445 
2s ae 6 O.LI4138 QO. 355084 
26> 715 0.260326 0.088162 
0.073026 0.435456] 
0.383268 O.3 71400 
0.416028 Ooh? 5 
0.225861 0.320409 
0.284265 
U.333147 
0.226094 
0.314591 
U.35304]1 
QO, 2 31555 
Osh 735 
0.285355 
0.201 147 
0.326326 

SR OHR 7HR 9HR 
0.165310 Oe Gilly 7 OS89G2 7 22 1994 
Ore) 0.271608 0.170280 Veo oT 

RANGE DISRWRIBCTION #2 

RADAR OG 3HR SHIR 
Oe e22 | Ora 9 N51 Orie. tS 0.347947 
Q, 244974 0.443105 0.312054 0.53/89 12 
0.220 2 Q.145520 0.255295 0.118954 
0.244750 O.141$40 0.339425 
0.195385 0.258664 0.084505 
0.063498 0.418551 
Dae 625 0.387630 
0.432849 0.142364 
0.231622 0.303028 
0.268500 
0.345476 
0.225661 
Ol a3@825 
Orso 33 
0.258948 
Ouro F | 
0.276574 
0.225085 
0.35428 7 

SHR 6HR 7HR O Islin 
0.150983 0.251506 0.332766 0.219481 
Onion? 0.296858 0.151890 0.285708 


TABLE 10 
TWO SAMPLE Tes neo cient 


RANGE DISTRIBUTION #1 


Comparison Reject 
( YeSa)) 
Radar vs 2hr balloon Ne oo 
Radar vs 3hr balloon NG 16 
Radar vs thr balloon NO 4 
Radar vs Shr balloon NO IO 
Radar vs 6hr balloon NO ol 
Radar vs 7hr balloon NO 49 
Radar vs 9hr balloon NO 2 
Radar vs > 4hr balloon No il 
RANGE DISTRIBU TiO s=2 
| Comparison Reject a“ 
(Yesi No) 

Radar vs 2hr balloon No aD 
Radar vs 3hr balloon Nie aS 
Radar vs thr balloon No 38 
Radar vs Shr balloon No 97 
Radar vs 6hr balloon NO 2 
Radar vs 7hr balloon Oo . 
Radar vs 9hr balloon NO sO2 
Radar vs > 4hr balloon 0 So 


Additionally, similar tests were conducted on data produced by appending 
current radar wind values to 2hr and 4hr stale balloon atmospheres (temperature and 
pressure). This data was then tested against the 2hr and 4hr stale balloon MET P(k)’s. 
These tests were to determine if the composite two or four hour data sets were as good 
as the two or four stale balloon data sets. In both cases, rejection of the null 
hypothesis that the composite data set was as good as the stale data set and 
acceptance of the alternative hypothesis was not possible, with exceptionally high a’ 


values required to reject. 
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Grobe at | 
teeeowealaeNeyY TEST RESELTS 


hee ae DISTRIBUTION 21 


Comparison Reject Ho ve 
(ES NO) 
Radar vs 2hr balloon NO 24355 
Radar vs 3hr balloon NO 0756 
Radar vs thr balloon Po 388 1 
Radar vs 3hr balloon NO 8467 
Radar vs 6hr balloon No yes. 
Radar vs /hr balloon NO 5370 
Radar vs 9hr balloon NO S000 
Radar vs > Shr balloon Oo .6837 
RaaGe DISTRIBUTION 22 
Comparison AReiect Ho o 
O) 
Radar vs 2hr balloon NO 3653 
Radar vs 3hr balloon NO 1463 
Radar vs thr balloon sO 4435 
Radar vs Shr balloon NO I184 
Radar vs 6hr balloon NO 4287 
Radar vs 7hr balloon io 6790 
Radar vs 9hr balloon NO 09 
Radar vs > 4hr balloon NO 8676 


Conover’s normal approximation, [Ref. 20: p. 217], was used for a’ > .5000 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONGEUSIONS 

The original data was very sparse and not gathered with this analvsis in nund. 
Hence, the results are not as conclusive as an experiment designed for that purpose 
could show. However, the objective of this thesis was to evaluate the results of a “pilot 
study” to determine if it was advisable to proceed with further testing and development 
of the WINDFINDER capabilities. The results of this analysis do in fact strongly 
support further testing and development. The large level of significance necessary to 
reject the null hypothesis in all comparisons for both distributions considered shows 
that not onlv can it not be concluded that the radar MET P({k) is worse than the stale 
balloon MET P(k), but there appears to be some evidence that the radar MET P(k) is 
as good as the stale balloon MET P(k). If this capability were available to the artillery 
units in the field during combat, it 1s conceivable that there would be a significant 
increase in the accuracy of first round artillery shots under rapidly changing weather 
conditions. 

This capability does not come for free. By using the FIREFINDER radar in the 
WINDFINDER mode, an increase in the electromagnetic radiating time 1s incurred. 
This increase in radiating time would most hkely result in some increase in the radar 
vulnerability, which could degrade the performance of the radar in its counter-battery 
mission. This issue of increased vulnerability versus increased first round accuracy 
during severe climatological conditions is of paramount concern to the artillery 
community as a whole. However, it is beyond the scope of this thesis, and is 


mentioned here only in the interest of completeness. 


B. RECOMMENDATIONS 

Given such strong indications of parity between the two MET svstems, the next 
and most logical, phase of testing should be pursued. These results warrant the 
expense of a specially designed and administered field experiment where both MET 
svstems are employed and targets are engaged by artillery of all calibers under as varied 
conditions as possible. This experiment should include normal, adverse, and rapidly 
changing Weather conditions, and be conducted during all times of the day and night, 
for it 1s then that the true strengths and weaknesses of both systems could be 


evaluated. 
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APE ENDIN A 
COMPUTER MET MESSAGES 


This appendix contains two computer MET messages. They are the data 
gathered by ASL on 28 February 1986 at O500hrs. The first message is the balloon 
MIET data, the second 1s the radar MET data. 

The columns are: 

e Altitude zone - these are explained in FM 6-15, [Ref. 4] 

e Direction the wind is conung from in mils x .1, i.e. 710 muls is listed as 71.0 
e Wind speed in knots 

@) Pemperatiive in t& «10:0, 1.c. 290.3 is hsted as 2903.0 

e Barometric Pressure in nullibars 

Since the radar is not effective in sensing the wind signature above the planetarv 
boundary laver (PBL), the wind values above the PBL must be estimated. The method 
used to estimate these values is to assume that they remain constant above the highest 
zone that the radar can sense, and set them equal to the values in the highest zone 
sensed. The analvtic atmosphere estimates temperature and pressure values within 
each zone, and 1s not limited by the PBL. 

The dates and ttmes that MET data was gathered at Yuma are listed below. 

poe nla slOS6 4 March 1986 


0500hrs O500hrs 
OSOOhrs OSOOhrs 
1400hrs 1000hrs 
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APPENDIN B 
INDIVIDUAL RANGE P(K) RESULTS 


This appendix is a listing of the calculated P(k)’s for the individual ranges. The 
first column is for 1Okm,. the second for 12km. and the third column contains the 
results for 1¢4km. The data are divided bv the tvpe of comparison indicated. 
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